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ABSTRACT 

An  asymptotic  method  is  developed  for  the  solution  of  linearised 
equations  governing  the  flow  of  a  compressible  viscous  fluid  with  a  free 
surface.  The  approach  used  here  is  based  upon  the  ray  method  expansion 
originally  developed  by  Keller.  A  general  uniform  asymptotic  expansion  is 
also  constructed  to  remove  anomalies  in  which  an  amplitude  function  in  the  ray 
method  expansion  becomes  infinite. 
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SIGNIFICANCE  AND  EXPLANATION 


-J  r,  4.  r  /  ' 

We  "develop  an  asymptotic  method  to  solve  the  linearised  Navi er- Stokes 
equations  governing  the  flow  of  a  compressible  viscous  fluid  subject  to  free 
surface  and  rigid  bottom  boundary  conditions.  The  solution  of  these  equations 
is  assumed  to  consist  of  a  phase  function  and  an  amplitude  function.  It  is 
found  that  the  phase  function  satisfies  the  Hamilton- Jacobi  equation,  and  the 
first  order  approximation  to  the  amplitude  function  satisfies  a  transport 
equation.  The  Hamilton- Jacobi  equation  may  be  solved  by  means  of  the  method 
of  characteristics,  which  reduces  the  equation  to  a  set  of  ordinary 
differential  equations.  Their  solutions  determine  a  family  of  time-space 
curves  called  rays.  Hie  transport  equation  can  be  easily  integrated  along 
each  ray  to  yield  the  so-called  conservation  relation.  At  certain  anomalies 
the  amplitude  function  becomes  infinite  and  a  uniform  expansion  is  then 
constructed  to  remove  these  difficulties. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


RAY  METHOD  FOR  FLOW  OF  A  COMPRESSIBLE  VISCOUS  FLUID 

M.  C.  Shan* 

I.  Introduction 

Since  the  ray  method  was  systematically  developed  by  Keller  and  his  co- 
workers,  much  progress  has  been  made  in  recent  years  and  a  review  of  the 
method  can  be  found  in  Keller  [1].  Although  the  ray  method  was  first 
formulated  for  electromagnetic  waves  and  has  been  extended  to  many  other 
related  areas,  it  has  also  found  a  lot  of  interesting  and  fruitful 
applications  to  the  wave  propagation  problems  of  a  fluid,  especially  with  a 
free  surface.  In  the  past  decade,  many  results  have  been  obtained  for  the 
case  of  an  irrotational,  inviscid  or  incompressible,  viscous  fluid  with  free 
surface.  However,  only  recently  the  Navier-Stokas  equations  governing  the 
motion  of  a  compressible  viscous  fluid  have  received  much  attention,  as  shown 
in  the  review  article  by  Solonnikov  and  Kazhiklov  [2] .  At  present  very  few 
solutions  of  the  linearized  Navier-Stokes  equations  subject  to  free  surface 
conditions  are  available,  to  say  nothing  of  the  nonlinear  ones.  It  should  be 
of  great  interest  to  develop  an  asymptotic  method  to  study  these  equations. 
Indeed  their  physical  applications  abound,  for  example,  wave  propagation  in 
the  atmosphere  or  ocean.  A  compressible  viscous  fluid  may  be  also  used  to 
model  flow  in  porous  media,  a  discussion  of  which  may  be  found  in  Bear  [3]. 

In  the  derivation  of  the  ray  method  by  Keller  [1]  for  surface  waves  on 
water  over  a  variable  bottom,  upon  which  our  approach  is  based,  a  wave-like 
solution  with  an  amplitude  function  and  a  phase  function  was  assumed,  and  the 
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boundary  conditions  at  the  free  surface  and  the  rigid  bottom  were  taken  into 
account  as  an  integral  part  of  the  formulation  of  the  method.  The  three- 
dimensional  problem  is  reduced  to  the  problem  of  constructing  rays  and  wave 
fronts  on  the  two-dimensional  equilibrium  free  surface,  and  the  solution  of 
partial  differential  equations  is  reduced  to  the  Integration  of  ordinary 
differential  equations  along  rays.  This  method  has  since  proven  very 
successful  in  dealing  with  the  problems  of  wave  propagation  in  a  fluid  with  a 
free  surface  (Shen,  [5]).  For  the  problem  considered  here,  we  adopt  the  form 
for  the  ray  method  expansion  used  in  Cohen  and  Lewis  [6] ,  Voronka  and  Keller 
[7] .  It  is  also  a  well  known  fact  that  the  ray  method  expansion  in  its  usual 
form  may  fail  at  anomalies  where  the  amplitude  function  becomes  infinite.  We 
then  use  a  modified  version  of  the  uniform  asymptotic  expansion  in  Shen  and 
Keller  [8]  to  remove  these  difficulties.  The  uniform  asyptotic  expansion  at  a 
caustic  was  first  developed  by  Ludwig  [9]  and  Krautsov  [10]. 

The  approach  used  here  may  be  sketched  as  follows:  We  consider  a  shallow 
layer  of  a  compressible  viscous  fluid  supported  below  by  a  rigid  bottom  and 
bounded  above  by  a  free  surface.  The  rigid  bottom  is  assumed  to  vary  slowly 
in  two  horisontal  directions.  A  large  parameter  (5  is  introduced  as  the 
ratio  of  the  horisontal  length  scale  L  to  the  vertical  length  scale  H.  The 
governing  equations  are  nondimensionalized  in  terms  of  the  dimensionless 
variables  appropriately  chosen.  Furthermore,  we  also  assume  that  the  constant 
viscosity  is  of  small  order  so  that  a  Reynolds  number  to  be  defined  is 
large.  First  we  construct  the  usual  ray  method  expansion  involving  a  phase 
function  and  an  amplitude  function.  By  substituting  the  expansion  in  the 
equations,  we  obtain  the  Hamilton-Jacobi  equation  for  the  phase  function, 
which  in  turn  determines  a  family  of  bi characteristics  called  rays.  The 
successive  approximations  to  the  amplitude  function  may  be  obtained  by 


integrating  transport  aquations  along  the  ray.  The  construction  of  a  uniform 
asymptotic  expansion  is  motivated  by  the  fora  used  in  [9],  [10],  and  based 
upon  an  ordinary  differential  equation  as  a  comparison  equation.  A  transfor¬ 
mation  is  discovered  to  reduce  the  results  for  the  ray  method  expansion  to 
those  for  the  uniform  asuaptotic  expansion.  By  appropriately  chosing  the 
coefficients  in  the  comparison  equation  and  a  conservation  relation  along  each 
ray  the  amplitude  function  becomes  finite  at  an  anomally. 

We  formulate  our  problem  and  construct  the  ray  method  expansion  in  $2. 

The  uniform  asymptotic  expansion  is  given  in  $3.  Discussions  about  some 
special  cases  and  possible  extension  of  the  results  obtained  are  presented  in 
§4.  Some  of  the  detailed  derivations  are  deferred  to  the  appendices. 
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2.  Hay  Method  Expansion 

The  linearized  Navier-Stokes  equations  governing  the  notion  of  a 
compressible  viscous  fluid  are  assumed  to  be  the  following: 

P#  *  +  <Pou#)  *  +  (P0V#)  *  +  (P0V)  *  “  0  ' 

t  x  y  z 

p.u  #  -  -p  *  +  P  +  (p/3)(V#  •  u  )  t  , 

0  t  x  x 

*  *  *  _2  *  — * 

p  v  #  *  -p  *  +  P  +  (p/3)(Vt  •  u  )  *  , 

t  y  y 

*  *  *  *  5  *  — * 

p.w  #  -  -p  *-pg+P  +  ( p/3 ) ( •  u  )  #  , 

t  z  z 

p*  -  (dp0/dpQ)p  *  (5) 

subject  to  the  boundary  conditions: 

At  z  -  0, 

*  * 

u  *  +  w  *  -  0  , 

z  z 


(1) 

(2) 

(3) 

(4) 


(6) 


*  * 

V  *  +  w  *  - 
z  y 


0  , 


(7) 


*  *  *  — * 

-p  +  p  gn  -  ( 2  p/3 )  V#  •  u  +  2yw  -  0 
u  z 


(8) 


n 


« 

t 


-  w 


0 


» 


*  *  •  *  — * 

at  z  *  -h  (x  ,y  ),  u  -  0 


(9) 

(10) 


Here  x*,y*  are  horizontal  coordinates  and  z *  is  positive  upward,  t*  is 

the  time,  V#  «  (3/3x*,  3^3y*'  3^3 u  "  (u*,v*,w#)  is  the  velocity,  pQ  is 

* 

the  equilibrium!  density,  considered  as  a  function  of  p^,  the  equilibrium 

*  * 

pressure,  p  is  the  density,  g  is  the  constant  gravitational  acceleration, 

#  * 

p  is  the  constant  viscosity,  z  *  n  is  the  equation  of  the  free  surface 
and  z*  ”  -h  is  the  equation  of  the  rigid  bottom  assumed  to  vary  slowly  in 
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« 


« 


*  *  •  « 
the  x  and  y  directions.  W*  not*  that  pp  t*  a  function  of  s 

•  *  *  * 
satisfying  P0  *  *  -Pq.  P0  "  0  •*  *  *  °* 

* 

Motivated  by  th*  assumptions  mad*  before,  we  introduce  nondimsnsional 
variables  and  parameters  as  below t 

(x,y,s)  -  (B~2x*,B”2y  ,a*)/H,  t  -  B-2t  /(H/g)^2  , 

q-  (u,v)  -  (u*,vV(gH)^2  ,  w  -  Bw*/(gH)^2  , 

•  *  *  *  *  * 

P  *  P  /lP0<0)gH],  pQ  -  p0/[p0(0)gH]  #  p  -  p  /P0(0)  , 

•  *  a  a 

P0  "  p0/p0(0)#  * 

a  1/, 

R  -  P0(0)(gH)/2H/M,  6  -  L/H  , 

where  we  aasuaa  R  is  large  and  without  lose  of  generality  we  set  R  *  3.  In 
terms  of  the  nondimenslonal  variables  (1)  to  (10)  becoete 

pt  +  V  *  <V>  +  0(POM*s  "  0  »  (11) 

p0^t  +  7P  *  B"3^q  +  ♦  (1/3)7(7»q)  +  (1/3)  ,  (12) 

Pn*  +  &3<P,+P>  *  8-3V2W  +  BW  _  +  ( 1/3)  ( 7»q)  +  ( 8/3 )  W  ,  (13) 

m  w  Jf  Zs  X  XX 

p  “  p0(p0)p  ’  (14) 


at  *  ■  0, 

B^q^  +  VW  ■  0  , 

B2(p-n)  +  (2/3) (fl”1 V  •  q  ♦  Wg)  -  2wz  -  0  , 

nt  -  BW  *  0  ; 


(15) 

(16) 
(17) 


at  z  m  -h. 


where 


q  ■  0,  w  ■  0  , 

7  -  <3/3x,  3/3y)  . 


(18) 


5 


For  ltrgt  0,  we  taiiM  that  tha  solutions  q,  W,  p,  p  and  n  possess 
an  expansion  of  the  fora 

4  -  exp [-3  S(t,x,y)](#1  +  6  %2  +  8  +•••)  ,  (19) 

Where  8  is  called  the  phase  function.  Substituting  (19)  in  ( ' 1  to  (18),  we 
obtain  a  sequence  of  equations  and  boundary  conditions  fox  succ;s.jive 
approxiaations  to  the  amplitude  function.  The  equations  for  the  irst 
approxiaation  are 


«Pl  -  p0k  .  q  +  (P0iya  -  o  , 

(20) 

-A  ‘  *pi  “  5iz.  ' 

(21) 

p1z  ■  "P1  ' 

(22) 

P1  "  PJ(p0)p1  1 

(23) 

at  s  ■  0, 

• 

at  z  -  -h. 

*1*  "  °'  p1  *  V  “ni  ’  W1  "  0  * 

(24) 

q,  -  0,  W1  -  0  , 

(25) 

where 

u  “  -S. ,  Jc  -  VS  . 
t 

From  (22)  to 

(24),  it  is  obtained  that 

p1 

-  n,  exp[-/J  pj(p0)dz)  -  Vo'  p,  -  pjp.n1  , 

(26) 

where  p^  -  -pQ, 

p  (0)  ”  1  have  been  used.  Now  let 

q  ■  , 

(27) 

and  it  follows  frca  (21),  (24)  and  (25)  that 

C1zz  ”  P0UC1  ”P0  ' 

(28) 

Q1b  -  0  at  z  -  0  , 

(29) 

Q1  -  0  at  z  *  -h  . 

(30) 

Assusie  that  the  Green's  function  of  (28)  to  (30)  is  given  by  G(z,z' ,u,h) . 
Then 
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(31) 


Cj  -  -/°h  PgGUjS'.w.Mds'  . 

Finally,  we  integrate  (20)  and  Bake  use  of  (25)  to  obtain 

-  “^(Pq  -  P0(-h>)  +  k\  flh  P0Q.,d*  ,  (32) 

where  k  ■  |k|.  However,  w1  must  also  satisfy  (24)  and  it  follows  that, 
assuaing  ?  0, 

wp0(-h)  -  k2  /®h  PgQjdz,  or  k2  *  H(u>,h)  ,  (33) 

which  yields  an  aquation  in  the  form  of  the  Hamilton- Jacobi  equation,  and  may 
be  solved  by  the  method  of  Characteristics.  The  corresponding  characteristic 
equations  are 


dt/do  ■  pH^,  dr/d o  -  2yk, 


du/do 


(34) 


dk/do  "  ”mVh'  ds/d0  “  P(2k  -  u>)  , 


where  u  is  a  proportionality  factor  and  r  ”  (x,y).  The  solutions  of  (34) 
determine  a  two-parameter  family  of  bicharacteristics 

(35) 

which  will  be  called  rays.  From  (34)  we  also  observe  that  to  is  constant 


t  -  t(o,o1,o2),  r  «  t(a,o^,a2)  , 


along  a  ray,  and  the  phase  function  is  determined  by 

S  -  S(o„)  +  f®  u(2k2  -  u)do  . 

0  ao 


(36) 


where  S(oQ )  is  the  value  of  S  at  some  point  o 


0g  on  a  ray. 


To  determine  we  proceed  to  the  equations  for  the  second 


approximation 


»P2  +  P,t  -  P0k  •  q2  +  +  (p0W2)s  -  0  , 


»P0q2  +  P0qlt  ‘  XP2  +  7pl  * 

P2z  *  _P2 


q2w  -  (1/3)W1b  . 


at  s  ■  0, 


2s 


p2  "  p£<Pq>P2  ' 


o,  p2  -  n2,  u»n2  +  n1t  -  h2  -  o 


(37) 

(38) 

(39) 

(40) 

(41) 


at 


*  -h, 


0 


(42) 


*2-°'  *2- 
Froa  (39)  to  (41),  we  obta  n  as  before  that 

p.  -  P0h2,  P2  -  Pfl<P0)P0n2 


Now  we  let 


i2  -  (kn2  -  +  Q2  , 


(43) 


(44) 


where  Q.,  satisfies  (28)  to  (30).  Then  it  follows  frost  (38),  (41)  and  (42) 
that 

®2«  '  “P0®2  “  p0<1t  +  {1/3)kW1z  '  (45) 


Q2z  -  0,  at  z  -  0 
Q2  *  0,  at  z  *  -h 


(46) 

(47) 


Integrating  (37)  with  respect  to  z  from  z  -  -h  to  z  «  0  and  making  use 
of  (27),  (33),  (41)  and  (42),  we  have 


/-h  [nltp0(p0,p0  +  P0Q1*  *  Vni  “  p0*  *  *2  +  V,(p0kVl)ldz 


(48) 


+  nit  "  °  ‘ 


Upon  multiplying  (48)  by  2n^  and  rearranging  the  terms,  we  finally  obtain 

(In2  t  +  7*(Ini  ~  An*  -  0  ,  (49) 

where 

I  -  P„(-h)  -  k2  /®h  ppQ,!*)  dz  , 

A  -  k2  J®h  pQ(z)  /®h  G(z,z\«,h)  [«op0z'  +  k2Q1(z')P0(z,)]dz*  . 

The  detailed  drivation  of  (49)  is  deferred  to  Appendix  A.  Let  J(o)  be  the 
Jacobian  of  transformation  from  the  ray  coordinates  o,  ,  c>2  to  the  t,  x, 
y  coordinates,  where  we  choose  t  -  a.  By  the  well-known  identity 


(49)  reduces  to 


J-1  dJ/dt  -  V» (dr/dt)  , 


d(IJn2 )/dt  -  AJn2  . 


(50) 


(51) 


Integration  of  (50)  along  a  ray  yields 

Jin?  *j cpI-J?  AI_1dt)  -  constant  . 

*0 

Aa  Man  froa  (51)  may  became  infinite  at  anomalies  where,  for  example, 

j(t)  -  0.  To  raaove  thee,  we  ahall  construct  a  uniform  asymptotic  expansion 
in  the  next  section. 


I 


3.  Uniform  asymptotic  expansion 


lhe  construction  of  a  uni fora  asymptotic  expansion  is  based  upon  the 
approach  developed  in  [11]/  and  a  discussion  of  the  aethod  aay  be  found  in 
[8] .  Assuae  that  q,  w,  p,  p  and  n  possess  an  asymptotic  expansion  of  the 
fora 

+  -  exp(-0e)[*(1>V<6VC)  +  3V“%(2)v'(8V5>]  ,  (52) 

v 

where  V(B  C)  satisfies  the  second  order  ordinary  differential  equation 


v-(8vC)  +  8"VP(5)v*(8VC)  -  B2"2V(5)V(8ve>  -  0  , 

♦  (1>  ~  ♦»>  ♦  ♦  B‘2^>  /  i  -  1.2  . 


(53) 


0  and  £  are  functions  of  t  and  r  only,  v  is  a  nonnegative  number,  and 

P(€),  Q(S)  are  functions  to  be  chosen.  We  define 

♦*  -  ♦(1,±Q*(2>  -  ♦*  ♦  B'%*  +  S'2**  ♦•••  ,  (54) 

8*  -  «  t  ll  Q(C')d*',  0*  -  ±(Q*  +  QP)  ,  (55) 

and  substitute  (52)  for  q,  w,  p,  p  and  T|  in  (10)  to  (17).  Next  we  collect 

V  V 

the  coefficients  of  V(B  €)  and  V'(B  C)  and  set  them  to  zero.  Then  we 


multiply  the  coefficients  of  V*  by  QB 


-v*1 


and  add  them  together  and 


subtract  them  from  the  corresponding  coefficients  of  V.  Making  use  of  (54) 


and  (55),  and  emitting  the  superscripts 


we  obtain 


B(pu>  -  pQq*k  +  w^)  +  Pt  ♦  V»(p0q)  -  D(p(2,?t  +  pQq(2)  »V£)  *  0 

"  «*>  +  poS:  +  7p  “  D(p0q(2>5t  +  P<2>75) 

-  Bq._  -  ( 1/3)kw  +  0<B_1)  . 

zz  z 

P3  +  p  •  o(B-1)  , 

p±  “  po(po)p  1 


(56) 

(57) 

(58) 

(59) 


at  a  •  0, 


(60) 


q^  +  0(6  2)  -  0,  p  -  n  +  0(6  2)  ■  o  , 

(2) 

6(«n  -  w)  ♦  n  -  D5to  «o  , 


at  z  ■  -h, 

q*  -  0,  w*  -  0  ,  (61) 

.1-2  -1  -2 
where  0(6  ) ,  0(6  )  denote  tens  of  orders  6  and  6  respectively  and 

will  not  be  used  in  the  sequel. 

The  equations  for  the  first  approximation  are 


“pi  ‘  pok  *  qi  +  tpoV*  "  0  ' 

(62) 

V*.  "  kp1  -  -1Z,  ' 

(63) 

P1*  "  -p1  ' 

(64) 

P1  "  p0(p0)p1  » 

(65) 

at  z  «  0, 

ii*  -  o,  Pl  -  iy  -  wt  -  o  , 

(66) 

at  z  -  -h. 

qt  -  0,  w1  -  0  . 

(67) 

which  are  the  same  as  (19)  to  (24).  The  results  obtained  for  the  first 
approximation  in  the  ray  method  expansion  can  be  carried  over  without  change, 
and  the  Kami  lton-Jacobl  equations  for  S*  are 

w*P0(-h)  -  (k*)2  J°h  pQdz  ,  (68) 

the  solutions  of  which  determine  two  two-parameter  families  of  rays 

t*  »  t±(o,o1,o2),  r*  -  r*(o,alfa2)  .  (69) 

To  obtain  a  conservation  relation  for  we  need  the  equations  for  the 

second  approximation 

—  —  —  (21 
“P2  +  P1t  "  PoVk  +  7*(poV  +  (p0w2’«  "  D(P1  et 

(70) 

+  P0q(2>‘?e)  -  o  . 

P0«52  ♦  p.q^  -  kp2  +  7pt  -  O(P0qS2Ut  +  p52,VO 

(71) 

•q2„-  (1/3)kl»u  . 
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“j.  *  «j  ’  0  • 

P,  -  PjlPo’.J  ’ 

(72) 

at  *  ■  0, 

« 2 .  ■  °*  ”2  •  n2  ■ 

(73) 

(2) 

“n2  -  *2  +  nit  ’  D5tni  "  °  * 

(74) 

a> 

ft 

N 

1 

ir 

% 

q2  “  °'  W2  "  °  * 

(75) 

From  (72) 

to  (73),  we  easily  obtain  that 

P2  "  poV  P2  "  p0(p0)p0n2  * 

(76) 

As  seen  from  (71),  (73)  and  (75),  we  may  let 

«2  -  (Sn,  -  Vn,)e,  ♦  Qj  *  Bj 

(77) 

where  Q1 ,  C2 

satisfy  (27)  to  (29),  and  (45)  to  (47)  respectively,  and 

*3 

satisfies 

«3,.  -  wp0®3  *  -“'"o'9!2’'!  *  > 

(78) 

03,  -»  «  .  -  0  , 

(79) 

C3  -  0  at  s  *  -h  . 

(80) 

We  now  integrate  (70)  with  respect  to  *  from  s  •  -h  to  s  •  0,  and  sake 
use  of  (27),  (68),  (74)  and  (75)  to  obtain 


f-h  tnitp0(p0)p0  +  p0fi1*‘7ni  •  P0*’°2  +  7,{p0*niQ1)]d‘  +  nit 

(81) 

■  f-h  v-V'  ♦  /!„  d<»12\  *  * o«52,,,',d'  *  °V'!21  • 

Note  that  the  left  hand  side  of  (81)  is  the  sans  as  (48).  We  Multiply  both 
sides  of  (81)  by  2t\j»  and  follow  the  same  derivation  given  in  Appendix  A 
to  obtain 

<in^)t  ♦  7»2k  J°h  njp^ds  -  An2 

(82) 

-  /!„  X'Sot'W*  +  /!h  D(pJ2,5t  ♦  P0qS2)»7C)ds  ♦  2n1DCtnj2)  * 

Zn  Appendix  2,  we  show  that  if 
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n1  -  *«K1#  and  2R*  -  RDQ*1  ,  (83) 

Than  C 1  satisfies  the  mm  conservation  relation  In  the  ray  Method 

expansion,  and  It  follow*  that,  along  a  ray 

J±I*R”2(n*)2eicp[-/^  A*Cl2)“1dtl  •  conatant  .  (84) 

0 

Proa  (55)  and  (79),  we  find  that 

R  -  J^expKI/2)  /SlCMdC'L  P  -  tlog  (R2/Q)l*  •  (85) 

At  an  anoawly,  we  May  ehooaa  R  ao  that  R  tend  to  a  finite  Unit 

there.  P,  Q  and  u  should  be  suitably  chosen  to  Make  (53)  as  ainple  aa 
possible . 
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4.  Discussion 


To  apply  th«  prsssnt  results  to  flow  of  a  compressible  viscous  fluid  with 
a  fixed  upper  boundary,  we  simply  replace  the  boundary  conditions  at  z  *  0 
by  q  ■  0,  «  *  0,  at  z  -  H,  the  equation  of  the  upper  boundary,  and 
by  p1(H),  the  value  of  p1  at  z  -  H.  The  Hamilton- Jacobi  equation  assuees 
the  form 

«(P0(-h)  -  Pg(H) )  -  k2  /®h  Q^d*  , 

where  Q1  satisfies  (28)  with  *  0  at  z  m  -h,  H.  The  conservation 
relation  along  a  ray  is  still  the  sane  as  (50)  except  X  and  A  are  now 
given  by 

i  -  P0(-h)  -  p0(h)  -  k2  j|!h  Q1u)p0a*  , 

A  -  k2  j"h  P0(x)  /®h  G(*,*',«,h,H)(li>p0B,  +  k2Q^ (x*  )pQ (a*  )]dx*  . 

The  uniform  asymptotic  expansion  for  this  case  can  also  be  dealt  with 
accordingly. 

As  an  illustration  of  the  construction  of  a  uniform  asymptotic  expansion, 
we  consider  a  strictly  convex  caustic  where  J(t)  “  0.  Following  [9],  we 
choose  u  *  2/3,  Q  “  ,  and  R  »  From  (85),  P  -  0  and  (53)  becomes 

the  Airy  equation 

V*(82/3C)  -  B2/3C  V(B2/3C)  -  0  . 

If  there  are  two  caustics  corresponding  to  C  ■  0,  C  *  then  (53)  becomes 

the  equation  for  Weber  functions 

V(5>  -  e25(e1-ov(C)  -  0  , 

where  v  •  0,  R  ■  C  ^*(C  ^-C)  F  ■  0  and  W  “  R2.  The  contraction  of  a 
uniform  asymptotic  expansion  at  a  line  of  sero  depth  may  follow  the  procedure 
developed  in  [12],  and  we  shall  omit  the  derivation  here. 
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Finally  we  make  mm  r Marks  regarding  the  possible  extensions  and 

applications  of  the  ray  method.  If  the  fluid  region  is  channel-like,  that  is, 

the  rigid  bottoai  varies  slowing  in  one  direction  only,  say  the  x* 

direction.  In  this  case,  we  introduce  the  following  nondimensional  variables 

_2  *  *  • 

(x,y,s)  -  (0  x  ,y  ,z  )/H  , 

q  -  <v,w)  -  P(v*,w*)/(gH)1/2  , 

and  the  remaining  ones  are  the  sane  as  before.  Then  we  assume  S  -  s(t,x) 
and  the  ray  method  can  be  carried  out  without  much  difficulty. 

In  this  paper  we  only  consider  the  case  of  large  Reynolds  number  R. 

With  a  change  of  the  stretched  variables,  the  ray  method  can  also  be  extended 
to  the  case  of  finite  R,  but  the  derivations  will  be  much  more  involved.  He 
also  note  that  the  results  obtained  are  only  valid  for  positive  phase 
function  8  and  then  the  nonlinear  terms  omitted  will  have  a  negligible 
effect.  One  of  the  interesting  applications  of  the  method  to  physical 
problems  is  to  consider  a  point  source  on  the  free  surface,  for  example  to 
model  a  well  in  an  unconfined  aquifer.  Mathematically,  we  simply  add  a  Dirac 
distribution  on  the  right-hand  side  of  (17).  First  we  construct  the  rays  from 
the  source  based  upon  (33)  and  (34).  The  phase  and  amplitude  functions  are 
given  by  (36)  and  (51).  Then  the  initial  conditions  for  them  are  determined 
by  the  asymptotic  expansion  of  the  exact  solution  of  a  canonical  problem  based 
upon  a  system  of  boundry  layer  equations  near  the  origin.  The  procedure  to 
solve  a  similar  problem  may  be  found  in  (11]. 


Appendix  A 


In  (48),  we  first  note  that 

2ni  f-h  poci  **7ni  +  2  ni7*(po*niQi)dl 


-  2 v«/®h  tn2p0 (?,<!*  , 

where  ■  0  at  z  «  -h.  Next  fro*  (27),  (45)  to  (47), 

Q2  *  /®h  G(*,*,,«,h)(p0(*,)(hn1Q1)t  +  n/3»cw1g]d*  . 

We  collect  the  terms  with  time  derivatives  In  (48)  to  obtain 


2niV 


III 


po(po,pod* 


+  «  -  2^i  n 


P0(X) 


ll 


G(z,z* ,u,h)dz'dz 

-  (h2)fc  P„ (“h)  +  2  /®h  k*(kn1Q1)t  ( s' )Pq (z* )dz* 

-  K  tVo(-h>  ♦  /!„  k2«?f0dl'l  • 


Where  we  have  changed  the  order  of  Integration  and  made  use  of  the  symmetry 
property  of  the  Green's  function  in  the  second  integral.  Hence,  (48)  assumes 
the  form 


3_ 

3t 


{n2[p0(-h)  +  /®h  k2Q2P(Jdz'}  ♦  2V  • 


kn2P0Qid* 


(A.1) 


Finally  we  go  back  to  (33).  Assuming  u  as  a  function  of  k,  we  have,  by 
differentiating  both  sides  of  (33)  with  respect  to  k^  where  k  »  (k^,k2). 

am/akj  -  ^i/dt  -  2k±  /®h  PjjQ^zIPjjt-h)^2  /®h  Q^dz]"1  .  (A.2) 

we  show  that 

*  f-h  p0Q?dE  "  /-h  *  (A‘3) 

From  (28)  to  (30),  we  find  that  satisfies 

Q1(i>zz  “  Po“®1w  “  P0C1  ' 

®1«z  *  °  '  at  «  "  0  , 

Ct(u  "  0  '  «t  z  -  -h  . 
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Hence , 

Xt  follows  that 
rO 


Q1«  “  ^  01  p0fiid*' 


(A.4) 


/"h  P0(*>Oi*,(*,d*  "  /!h  p0<*>  J-h  «<*»*' »“**‘>P0<*,>O1<‘,>d*‘ 


by  Integration  by  parts,  symmetry  of  the  Green's  function  and  (31).  From 
(A. 1)  to  (A.3),  (49)  follows. 
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Appendix  B 

MB  substitute  n1  ■  R(?)C,  in  (82)  end  obtain 

*2I<lC*>t  ♦  27«k  /®h  PoQ^d*  "  As2]  -  -2RR'[IC^t 

♦  27? »k  /®h  C^P^ds]  ♦  R2  /^h  2?1tPflk«Q3ds  +  (B.t) 

♦  D(*S2)et  ♦  P0q7<2>*  V5)  +  *>v!2>ld*  "  ' 

12)  +  — 

where  -  (?1  -  Z^)/{2q).  It  we  can  show  that  the  RHS  of  (B.1)  is  sero 

under  (83),  then  C 1  satisfies  the  sane  transport  equation  as  (49),  and  (84) 
follows. 

From  (26),  (27)  and  (54), 

7?  -  (k  “  -  k  +)/(2Q),  Ct  "  (H+  -  hf)/(2Q)  , 

pj2>  -  P0(n*  -  n‘)/(2Q),  pj2>  -  Pop0tni  *  »>|>/<2Q)  * 

q7(2)  -  (k  *t)*Q*  -  k  ‘n"Q")/(2Q)  . 

We  substitute  them  in  the  RHS  of  (B.1)  and  consider  the  upper  sign  first.  The 
coefficient  of  (C+)2  is 

(u  -  u  )i  /2  ♦  (k  -  k  )  •  k  /_h  p^Q^dx 
which  can  easily  be  shown  to  be  sero  by  use  of  (31)  and  (A. 4).  The 
coefficient  of  C* c  is  found  as 

k  +  •  k”  { (u>+-«  )/®h  pQ(s)  /®h  pQ(s* )G*(s,s' ,u+,h)Q~(s' )dz'ds 

(B.2) 

+  /®h  PQ(s)  /®h  P0(s,)(G+(s,s,,«+,h)  -  G~(z,s*  ,o>**,h)]ds'ds}  , 

where  we  have  used  the  relations 

w±P0(-h)  -  (k*)  /®h  pgQ^ds  . 

To  show  (B.2)  is  sero,  we  apply  (28)  and  (30)  for  Q^t 

*  Azz  ■  “  ’po  '  (B*3) 

fil,  *  0  *  ■  - 0  •  <»•«> 
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Q*  -  0  at  s  -  -h  .  (B.5) 

Than  1st  UVdx  “  (D,V) ,  and  we  have  from  (3) 

<Q>+q;>  -  +  (Q“,-P0<d+q'J) 

(B.6) 

-  (qV-p0«‘q;)  -  <Qv-P0>  -  <Q^-P0)  • 


By  tha  boundary  conditions  (B.3)  and  (B.4),  (b.5)  becomes 

(«+  -  »-)(p0Q*,q’>  +  (P0,Q^)  -  (P0>Q^)  -  o  . 
the  VB8  of  tha  above  aquation  reduces  to  (B.2)  by  integration  by  parts, 
symmetry  of  tha  Green' a  function  and  (31).  Similarly  we  can  also  show  that 
the  RHS  of  (B. 1)  coresponding  to  the  lower  sign  is  also  equal  to  sero. 
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